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Abstract 

The observation that Type la supernovae (SNe la) are fainter than expected given their red shifts 
has led to the conclusion that the expansion of the universe is accelerating. The widely accepted 
hypothesis is that this acceleration is caused by a cosmological constant or, more generally, some 
dark energy field that pervades the universe. In this paper, we explore what, on their own, the 
supernovae data tell us about this hypothesis. We do so by answering the following question: can 
these data be explained with a model in which the strength of gravity varies on a cosmic timescale? 
We conclude that they can. Consequently, the supernovae data alone are insufficient to distinguish 
between a model with a cosmological constant and one in which G varies. However, the varying-G 
models prove not to be viable when other data are taken into account. This topic is an ideal 
one for investigation by an undergraduate physics major because the entire chain of reasoning 
from models to data analysis is well within the mathematical and conceptual sophistication of a 
motivated undergraduate. 



I. INTRODUCTION 



Physical cosmology^ is the branch of astronomy concerned with the large scale struc- 



ture and evolution of the universe. Cosmology is unique in that two apparently disjoint 
disciplines — the physics of the very small and the physics of the very large — are both needed 
to achieve our current understanding of the universe. It is remarkable, for example, that 
quantum fluctuations that occurred on microscopic scales in the very early universe may 
have left an imprint on the largest structures in the universe. The observation that the 
universe contains matter — with only trace amounts of antimatter, rather than matter and 
antimatter in equal amounts, may find its explanation, one day, using earthbound particle 
accelerators. Dark matter, for which there is much compelling evidence,^ may yet turn out 
to comprise weakly interacting particles that may be accessible in laboratories. The rela- 
tively recent synergy between the theories of the very small and the very large is a thrilling 
achievement. However, there is a cloud on the horizon called dark energypEl 

A big surprise in cosmology came in 1998 when the High-Z TeanP and the Supernova 
Cosmology ProjectP both observed, independently, that Type la supernovae were fainter 
than expected. After careful consideration of alternative explanations, both teams of re- 
searchers interpreted their observations as evidence that the SNe la are further away than 
expected, given their red shifts and assuming a decelerating universal expansion. If the SNe 
la are further away than expected, then the average expansion rate of the universe since the 
big bang must be higher than previously thought. Both teams, in fact, went further: they 
concluded that the expansion of the universe is accelerating. Today, the broadly accepted 
hypothesis is that this acceleration is driven by a form of energy called dark energy that 
pervades the universe. In the simplest model, dark energy is identified with the cosmological 
constant A that appears in the general form of Einstein's theory of gravity, general relativity. 
In more complicated models,- dark energy is modeled as a dynamical, evolving, field. 

Cosmologists have created a compelling and coherent cosmology based on the Friedmann 
equation 




(1) 



and the associated Friedmann-Lemaitre- Robertson- Walker (FLRW) metric^ , 



ds 2 = c 2 dt 2 - a 2 (t) 



dr 2 



+ r 2 {d0 2 + sin 2 dcj) 2 ) , 
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where a(t) is the dimensionless scale factor — normalized so that a (to) = 1 & t the present 
time t — to, a = da/dt, G is the gravitational constant, p is the density of all forms of 
energy^ excluding the contribution from the cosmological constant, A, and — oo < K < oo 
is the spatial curvature. The radial coordinate r is defined so that the proper area of a 
sphere, centered at any conveniently chosen origin, is Aq = 4irr 2 at the present time. As 
usual, symbols with a subscript of zero denote quantities evaluated at t — to- The comoving 
distance \ associated with the radial coordinate r is given by 

dr' 



X 



I 



y/l - Kr' 2 ' 

sm-^KWryK 1 ' 2 , (3) 



while d(t) = a(t)x is the proper distance at time t. By construction, the comoving and 
proper distances are numerically identical today. The radial coordinate r, comoving distance 
X, radius of curvature K' 1 ^ 2 , and proper distance d(i) are conventionally measured in Mega- 
parsecs (Mpc). Inverting Eq. (|3]), we obtain 

r = sm(K^ 2 X )/K^ 2 . (4) 

For a spatially flat universe, that is, one with K = 0, Eq. Q simplifies to r — X- 

The standard model of cosmology, with K = and A > 0, works remarkably well; 
however, current physics predicts^ a value of the cosmological constant A that exceeds the 
observed value by a factor of at least 10 50 ! This difficulty motivates the exploration of 
alternative explanations, such as ones that invoke time-varying "constants" P After all, we 
know of no compelling reasons why the parameters that appear in our current theories of 
the physical universe should be independent of space and time. From some perspectives, 
the puzzle is why they should be constant at allP 

Another motivation for exploring alternative explanations of the supernovae data is to de- 
termine whether they alone are sufficient to distinguish between a model with a cosmological 
constant and models without, such as the varying-G models we consider in this paper. 

A third, rather different motivation, is the pedagogical value of such investigations. This 
topic is ideally suited for directed individual study (DIS) by an undergraduate physics ma- 
jor. It is exciting and lends itself to open-ended exploration. The work reported below 
was undertaken by one of the authors (RD), an undergraduate physics major, under the 
supervision of the other (HBP). We share the view of many that, ideally, all undergraduates 



should be afforded the opportunity to engage in authentic research. However, this is not 
always easy: many exciting topics unfortunately require rather more material than can be 
mastered in a reasonable amount of time by a busy student. The advantage of cosmology 
is that it is intrinsically interesting to many students and, provided one chooses the topic 
carefully and one is prepared to make appropriate conceptual approximations, one can find 
interesting cosmological studies that can be done using mathematics and concepts that are 
readily accessible to an undergraduate student. We fully endorse the idea that, for such 
students, a "mathematics first" approach followed by applications is less desirable than the 
"physics first" approach as advocated by Hartle for general relativityP The cosmological 
investigation described below was done in that spirit. 

This paper explores two simple phenomenological models of varying-G cosmology 8 using 
the data compiled by Kowalski et al. 10 on 307 supernovae. We assume a spatially flat 
(K = 0) universe (motivated in part by the expectations from inflation^) and we set A = 0. 
However, for completeness, we write all expressions in a form that is valid for arbitrary 
values of K and A. 

We find fits to the supernovae data that are competitive with the simplest dark energy 
model. The fact that non-dark energy models can account for these data is a reminder 
that the supernovae data alone are insufficient to establish dark energy as the preferred 
hypothesis. That hypothesis becomes compelling only when different data-sets are analyzed 
together. Likewise, any varying-G model must fit not only the supernovae data, but must 
also be in accord with other data. However, in view of our expressed goal to provide 
an example of an authentic research project that can be conducted in its entirety by an 
undergraduate student, we restrict the scope to only one other datum: the bounds on G/G 
at our current epoch. We find that our two varying-G models fail the bounds on G/G, 
thereby ruling out this form of variation in G. An interesting aspect of the first model is 
that the scale factor becomes infinite in a finite amount of time. In such a model, all comes 
to an end in a catastrophic shredding of everything, a doomsday scenario that has been 
dubbed the big ripJ^ 

The paper is organized as follows. Section [IT] describes those additional aspects of FLRW 



cosmology that are needed to understand the SNe la data. Section HI motivates the varying- 
G Friedmann equation we have used to describe the evolution of the scale factor, while 



Sec. IV introduces our varying-G models and summarizes the results we obtained from 
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them. The paper ends with a discussion and concluding remarks. 



II. SUPERNOVA COSMOLOGY 

A key problem in observational cosmology is measuring distances to galaxies. To do that 
we need two things: 1) a standard candle and 2) an operational definition of distance. We 
consider first the standard candle. 

A standard candle is a source whose absolute luminosity is known. Type la supernovae^ 
are currently the best "standardizable" candles for very large distances. A Type la supernova 
is believed to occur when a star in a binary system overflows its Roche lobe (the region 
within which its matter is gravitationally bound) causing material from it to accrete onto 
the companion white dwarf. The mass of the white dwarf gradually increases towards 
the Chandrasekhar limit (of about 1.4 solar masses),^ triggering runaway nuclear burning 
within the star that releases more energy in a matter of weeks than the sun will emit in 
ten billion years. In another class of models, an explosion is triggered by the merger of 
two low-mass white dwarfs. In a third class of models, a carbon-oxygen low-mass white 
dwarf explodes when the helium, accreted from a companion star, detonates. For a good 
review of Type la supernovae models, see RefP. By measuring specific characteristics of the 
supernovae light curves (graphs of brightness as a function of time), it is possible to make 
empirically-derived corrections for the observed variations in SNe la brightness and thereby 
create well-calibrated standard candlesP 



A. Luminosity Distance 

The proper distance between two points in space is a well-defined concept, but, unfortu- 
nately, it cannot be measured in practice. Instead, astronomers use a definition of distance 
based on the flux of energy received on earth from the luminous object, that is, on the 
energy received per unit area per unit time, 



where L = dE/dt is the object's luminosity, that is, its rate of total energy emission, and 
A = 4nr 2 is the proper area, at t — to, of the sphere centered at the location once occupied 
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by the supernova. This formula for the flux is valid for a static universe and for a source that 
emits energy isotropically. In an expanding universe, however, the luminosity L crossing this 
sphere is diminished by the factor (1 + z) 2 . One factor of 1 + z arises from the reduction in 
energy of each photon received on earth relative to the energy it had at emission, yielding 
dE — >■ dE/(l + z). By definition, the red shift z = (A r — A e )/A e where A e and A r are the 
emitted and received wavelengths, respectively. The second factor of 1 + z is due to the 
reduction in the rate of arrival of photons at the earth, which yields 1/dt — > (l/dt)/(l + z). 
The corrected expression for the flux is then 

f = 4tt[(1 + z)r] 2 = 4^f' ^ 

where dj, = (1 + z)r is called the luminosity distance. For arbitrary values of the curvature 
K, the radial coordinate r is related to the comoving distance \ v i & Eq. Q, which, as noted 
above, reduces to r = \ when K = 0. 



B. Distance Modulus 



Astronomers measure energy fluxes. But, by convention, fluxes are converted into magni- 
tudes m defined by / = ql0~ 2m ^ 5 = Lj (AndD, where q is the flux from objects of magnitude 
zero, with the luminosity distance di measured in Mega-parsecs (Mpc), and absolute mag- 
nitudes M defined by f M = glO" 2M/5 = L/(4nd 2 M ) where d M = 10" 5 Mpc. Next they take 
the logarithm to base 10 of the ratio f M /f = lO°- 4 ( m - A/ ) = (c^/KT 5 ) 2 to arrive at 

m-M = 51og 10 [d L /10- 5 ], 

= 51og 10 [(l + z)r(z)/l(n 5 ]. (7) 

Note that in the ratio the constant q cancels. The difference \x = m— M between the apparent 
magnitude m of a source and its absolute magnitude M is called the distance modulus. The 
analysis of a supernova light curve results ultimately in two measured quantities: /x and z. 
The dataP^ used in our study are plotted in Fig. [lip'. The cosmology is all contained in the 
dependence of the radial distance r, or, equivalently, the comoving distance \i on the red 
shift z. The red shift, in turn, is related to the dimensionless scale factor, a(t), as follows 

a = l/(l+*). (8) 
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FIG. 1: Measurements of (fj,, z) for 307 Type la supernovae, compiled by Kowalski et al. 



Given a functional relationship between the comoving distance \ and red shift z, the 
distance modulus function in Eq. ^ can be fitted to the data in Fig. [I] to extract the 
parameters of the cosmological model. A formula for the comoving distance x can De 
deduced from the FLRW metric, Eq. ([2]), by noting that light in vacuum travels on null 
worldlines (for which ds — 0). Therefore, a light ray from a supernova at red shift z satisfies 
the relationship cdt = a(t)d\- Hence, 

X(z) = cf ^ (9) 

Jil+z)- 1 aa 

noting that the light ray was emitted when the scale factor was a = 1/(1 + z) and received 
when it assumes the value unity, today. 
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III. A VARYING- G FRIEDMANN EQUATION 



Our first assumption, which we have alluded to, is that the universe has zero spatial 
curvature. Our second assumption is that the Friedmann equation, Eq. Q, for a K = A = 
universe remains valid when G is allowed to vary with time. By valid we mean that the 
equation is a good approximation to some (unknown) exact equation describing the evolution 
of the scale factor in a universe in which G varies. This is an example of a conceptual 
approximation that renders the problem tractable for an undergraduate student. If we wish 
to remain strictly within the framework of general relativity, we should be cautious about 
replacing Eq. ([!]) with one in which G is a function of time because the Friedmann equation 
is derived from Einstein's equations 

Gfj, v + Ag^ = -^j— T^„, (10) 

which do not permit variations^ in G. The tensors G^ u and T^ v are the components of the 
Einstein and energy-momentum tensors, respectively, and g is the metric tensor. In order 
to allow for a possible variation of G, theories more general than Einstein's are needed, such 
as the scalar-tensor theories^^ in which gravity is assumed to couple to a scalar field 0, 
which — for weak constant coupling — yields the relationship G oc This, in turn, yields a 
modified Friedmann equation with a time-dependent G and additional terms of order G/G. 
If the latter terms are small enough, we obtain a Friedmann equation identical in form to 
the standard one, but with a time-dependent G. 

Writing G(t) = Gof(a), where Go is the current value of G and f(a) describes the 
assumed dependence of G on the scale factor a(t) and therefore cosmic time t, and using 
the definitions 

p c0 = 3H*/8nG , 
p A == Ac 2 /8nG , 
tyitf(a) = p(a)/pco, 

= Pa/PcO, 

n ee n M (i) + n A , (11) 

where p c o is the critical density now, Qm{(i) is the matter density parameter, and Ho is the 
Hubble constant — that is, the value of the Hubble parameter H(t) = a/a today, we may 
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write the modified Friedmann equation as 

(^J =H^[f(a)Q M (a) + (l-Q )a- 2 + Q A ], (12) 

noting that 

-Kc 2 = H 2 (l-n ). (13) 

With these definitions, we can write the expressions for the comoving distance x( z ) an d the 
universal time t(a) as follows, 

X ^ ^ H J(i +z) -i a 2 a/ ' f(a)Vt M (a) + (1 - fi )a~ 2 + fi A ' ^ ^ 

and 

i( a ) = -L f ^ (15) 

-no Jo x^f(x)Q M (x) + (1 - fi )^ -2 + 

The lifetime of the universe is given by to = t(l). 

Our third assumption is that the total mass-energy in the universe, whatever its nature, 

scales in the same way as matter; that is, we assume that Q\ = and CIm(o) = VLmq/^i 

where Qmo — ^m(I) denotes the value of the matter density parameter today. Since we 



also assume K = 0, Eqs. (11) and (13) show that, necessarily, = Qmo = 1- However, 
observations, interpreted within the context of the standard cosmology^ indicate that the 
matter density parameter Qmo ~ 0.3. The difference between Q = 1 and is presumed 
to be due to the cosmological constant or dark energy. If we wished to be consistent with 
this value of fiMo, while keeping A = 0, we need to use a model with K < 0. 

One of our goals, however, is to ascertain whether the SNe data, on their own, are 
sufficient to conclude that the A > model is preferred. To do so, we need merely exhibit 
another model that works as well. Here we consider varying-G models with K = A = and 
therefore Qmo = 1- Alternatively, one could consider A = 0, K ^ 0, models. It should be 
noted, however, that the curvature term (1 — Q )a~ 2 cannot accelerate the expansion. In 
a universe dominated by curvature, the Friedmann equation is a = constant, which implies 
zero acceleration. To obtain acceleration, one needs a term that dilutes less rapidly than the 
curvature term, which is the case for a cosmological constant or for the varying-G models 
described below. 
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IV. VARYING- G MODELS AND RESULTS 



In principle, a model for the variation of G should arise from some deep theoryP This, 
however, is far beyond the scope of this paper, which is to give an example of an interesting 
cosmological study than can be executed in its entirety by an undergraduate physics major, 
but that nonetheless yields interesting results. We proceed in a purely phenomenological 
manner. Our basic premise is that the supernovae are further away than expected because 
gravity was weaker in the past and, consequently, the universe decelerated less rapidly than 
would be the case were G constant and equal to its current value, Gq. 



A. Fits to Supernovae Data 

We studied several forms for the function f(a) in G(a) = Gof(a), but in this paper we 
report results for only two of them, each with a single adjustable, dimensionless, parameter, 
b. The first varying-G model we studied is defined by 

f(a) = e b{a - l \ model 1. 

In this model, there is no limit to how strong gravity can become. Another model studied 
is defined by 

f{a) =2/(1 + e - b(a_1) ), model 2, 

in which G is limited to twice its current value in the distant future. We have normalized 
both models so that G(a) assumes its current value, Gq, when a = 1. For K = models, 
the distance modulus, Eq. (|7j), may be written as 

fi(z, b,Q) = 5 log 10 [(l + z)H r(z)/c] + Q, (16) 

where the offset Q determines the vertical location of the modulus curved Note that 
Hqt{z)/c is dimensionless and independent of the Hubble constant. 



Evaluating Eqs. (14) and (15) for model 1, with K = (that is, fio = 1) and = 0, we 
find 



and 



x ( z ) = r (z) = —e b / 2 ^Jb[eii{^Jbj2) - erf( ^6(1 + z)"72)], (17) 

t(a) = ^[e b/2 {^/2T^eii{^/abj2) - 2 v ^e' ab/2 )]/b. (18) 
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This model exhibits a striking feature: the scale factor becomes infinite in a finite amount 



of time. We shall return to this point below. For model 2, the integrals in Eqs. (14) and 
(15) are evaluated numerically using the mid-point ruleP^ 



We fit Eq. (16) to the SNe data in Fig. [I] by minimizing the function 

307 

Y^lPn - n{z n ,b,Q)] 2 /a 2 n , 



X 



(19) 



n=l 



th 

with respect to the parameters b and Q, where 1 < n < 307 labels the n supernova at red 
shift z n and distance modulus fi n , measured with an uncertainty of ±o"„. The minimization 



of Eq. (19) is done using the program TMinuit, which is part of the ROOT data analysis 



package from CERN.^ For model 1, we get the result shown in Fig. [2j The fit gives the 
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FIG. 2: Residuals, /x n - /j,(z n ), of fit of model 1 to SNe la data (x 2 /ND = 313.0/305 = 1.03). 
value b = 2.09 ±0.08, from which we infer a lifetime of to = 15.1 ±0.3 (70 km s _1 Mpc _1 /iJo) 
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GyrP^The fact that the \ 2 P er degree of freedom (x 2 /ND) is 1.03 suggests that the modulus 
uncertainties are estimated correctly and that model 1 provides an excellent description of 
the data.^ A similarly good fit is found for model 2, as shown in Fig. [31 This fit yields b = 



3r 
2r 
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FIG. 3: Residuals, [i n - fi(z n ), of fit of model 2 to SNe la data {x 2 /ND = 316/305 = 1.04). 

3.27±0.11, with a x 2 /ND = 316/305 = 1.04. We find t = 16.2±0.4 (70 km s^Mpc" 1 /^) 
Gyr. For the simplest dark energy model, for which f(a) = 1 and VLm{o) = (1 — ^a)/& 3 with 
Vl A > 0, we find fi A = 0.71 ± .02 and t = 14.0 ± 0.3 (70 km s^Mpc-V-Ho) Gyr, consistent 
with the accepted results.® The \ 2 P er degree of freedom of the fit is 310/305 = 1.02. 

Since there is no compelling statistical basis to reject any of these models, we conclude 
that the supernovae data alone are insufficient to distinguish between them. However, these 
data when analyzed along with others® are consistent with a simple cosmology in which dark 
energy mimics a cosmological constant with the value Q\ ~ 0.7. The varying-G models 
should likewise be analyzed along with other data to see if a consistent picture emerges. 
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The fact that a K = A = model requires Qmo — 1> while the preferred value from galaxy 
and galaxy cluster measurements is Qmo — 0.3, is already an indication of a difficulty A 
systematic analysis of the relevant data, however, is a large task beyond the scope of this 
paper. Instead, we illustrate the importance of including other data by comparing the 
predicted fractional variation of G, G/G at the present epoch, with the available bounds. 



B. Bounds on the Variation of G 

The possible variation of G is usually characterized by the quantity G/G, which in terms 
of the logarithmic derivative of the function f(a) is given by 

G _ 1 dG _ d\nf 

G G da da ' 

where we have used the fact that Hq = a/ a = a at the present epoch. Figure [4] shows G/G 
as a function of the scale factor, for models 1 and 2. We see that at a = 1, G/G is equal 
to 1.5 x 10~ 10 y _1 and 1.15 x 1CT 10 y _1 , respectively. Unfortunately, these values for G/G 
are one to three orders of magnitude larger than the upper bounds that range from about 
1CT 10 y _1 to 10~ 13 y _1 depending on the method used to extract the bound 



V. DISCUSSION AND CONCLUSION 



Engaging undergraduate students in research can be an effective way to keep them ex- 
cited about science. Cosmology is particularly well suited in this regard because it is pos- 
sible to find topics that are both manageable in scope and scientifically interesting. We 
have presented an investigation of varying-G cosmological models that serve as examples of 
interesting research problems that are well matched to the mathematical sophistication of 
an undergraduate. 

The two phenomenological models we presented, in which the strength of gravity in- 
creases with cosmic time, provide excellent fits to the Type la supernovae data. We therefore 
conclude that the supernovae data alone cannot establish the dark energy hypothesis un- 
ambiguously. Other data are needed to render this hypothesis plausible. However, both our 
varying-G models fail to satisfy the bounds on G/G. Consequently, the particular variation 
of G described by these models is ruled out. In fact, one can make a stronger statement: all 
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FIG. 4: Fractional change in G per year for model 1 (solid curve) and model 2 (dashed curve) as 
a function of the dimensionless scale factor a(t). 

varying-G models that give rise to accelerated expansion, and that are based on the FLRW 
metric and the Friedmann equation, are ruled out by these bounds.^ Consider, for example, 
matter-dominated models, for which the Friedmann equation is H 2 ~ G/a 3 — > G ~ ad 2 . 
This yields G/G ~ H + 2d/d, from which we conclude that G/G <^ H. But a value of G/G 
of the order of H is inconsistent with the bounds on G/G, which are less than H by one 
to three orders of magnitude. 

The inability to distinguish between models (model degeneracy) is inherent in the Fried- 
mann equation because the latter is sensitive only to the total energy density of the universe 
and is agnostic with respect to how the energy density arises. This may be seen by writing 
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the Friedmann equation, Eq. ([I]) 



as 



) 



2 



8ttG 



a 



(Pm + Pk + Pa) 



(21) 



a 



3 



where pm, Pk, and p\ are the density contributions from matter, the curvature and the 
cosmological constant, respectively. The Hubble parameter a/ a is related to the sum of 
Pm + Pk + Pa not the individual components; or, equivalently, to the sum 



Therefore, it is possible to entertain different interpretations of the total energy density. For 
example, any model based on the Friedmann equation can be reinterpreted as one in which 
matter, perhaps of several different sorts, is either created, destroyed, or both, as the universe 
evolves. Consider, for example, the simple cosmological constant dark energy model, for 
which f(a) = 1, Qq — 1 and the total energy density is given by Q = (1 — Q\)/a 3 + Q\. 
This can be rewritten as Q = Q' M (a)/a 3 with Q' M (a) = 1 — Q\ + Q^a 3 - Since a~ 3 is the 
dilution factor for matter, the function Q' M (a) describes an increasing matter density in 
a comoving volume, which can be interpreted as the creation of matter as the universe 
expands! Alternatively, as done here, one can maintain the mass continuity equation, in 
which case matter is neither created nor destroyed and Q = Q^o/a 3 , but allow G to vary 
like f(a) ~ 1 + a 3 (l — Qmo)/^mo- Because of the invariance of the Friedmann equation 
with respect to such changes in interpretation, it is necessary to impose constraints on the 
cosmological parameters to remove the model degeneracy. Such constraints can come from 
other data, or other equations, or both. 

It seems odd, at first, that the strengthening of gravity with time leads not to the eventual 
gravitational collapse of the universe, but rather to its accelerating expansion. The reason 
for this is that every form of energy contributes to the geometry of spacetime. A model in 
which the strength of gravity changes with time is equivalent to another model in which the 
energy density changes in a specific way. If the energy density dilutes more rapidly than 
a~ 2 , then the expansion will slow down. If the strength of gravity increases such that in the 
equivalent, constant-G model, the energy dilutes more slowly than a~ 2 , the expansion will 
accelerate. In our varying-G models, the effective energy density increases with time. 

For model 1, the increasing strength of gravity leads to a startling prediction: a catas- 
trophic end to such a universe. This conclusion follows from the limit a — > oo of the lifetime 



Q = f(a)Q M (a) + (1 - n )a~ 2 + Q A . 
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expression, Eq. (18). We find that 



t(a ^00) = ^- exp(6/2) ^jr/b/b. (22) 
Ho 

According to this model, the universe has a finite lifetime of about 33 Gyr and will tear itself 
to pieces in its final moments! Such behavior has been dubbed the big rip and is a feature 
of cosmological models containing phantom energy:' 11 Within regions that are dominated by 
non-gravitational forces, the effect of a cosmological constant does not change with time and 
consequently the accelerating universal expansion will not disrupt already bound systems. 
By contrast, as the universe ages the effect of phantom energy increases in any finite volume 
of space. Eventually, this precipitates an escalating cascade of destruction at ever smaller 
scales until everything is torn asunder. We can only hope that phantom energy is just that: 
a phantom! 
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able to be measure the flux across all wavelengths. In practice, however, fluxes are measured in 
wavelength bands defined by standard filters, such as the 5-band filters.^ An observed super- 
nova spectrum is red shifted with respect to the spectrum in the rest frame of the supernova. 
Therefore, in general, a filter will transmit a flux that differs from that which would have been 
measured were it possible to measure the flux in the supernova's rest frame. Astronomers use 
corrections called K corrections to map the measured flux back to its value in the object's rest 
frame. Given a model of the object's spectrum, it is then possible, in principle, to infer the 
bolometric flux and hence the bolometric magnitude of the object. The distance moduli data 
compiled by Kowalski et al. are derived from l?-band magnitude data. 

The absolute luminosity of a supernova cannot be determined independently of the Hubble 
constant Hq. Consequently, in the fit of the modulus function to the data, it is only the shape of 
the function that contains useful information about the cosmology. The offset Q depends both 
on Hq as well as on the flux corrections. 

If the interval [a, b] is divided into N intervals of width h = {b — a)/N, the mid-point rule is 
£/(*)ds«fcl£ 1 /(a+(i-0.5)/0. 

The lifetimes can be computed given values for the parameters b and Hq. However, since we 
cannot extract a value of Hq from the fits, we compute the lifetimes using the nominal value 70 
km s _1 Mpc _1 for the Hubble constant, but we write all lifetimes in terms of the parameter Hq 
to make clear how the numerical values will change if Hq differs from the nominal value. 
If the reported modulus uncertainties are Gaussian distributed, one expects x 2 to be sampled 
from a probability density with mean N — P = ND, where ./V = 307 is the number of data 
points and P ~ 2 is the number of adjustable parameters. Therefore, for a fit that neither 
over-fits nor under-fits, we expect < x 2 > /ND ~ 1. The quantity P would be exactly equal to 
2 if the constraints that define the parameter estimates were linear in the parameters. 
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